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Abstract 

Waveguides of subwavelength dimensions with asymmetric geometries, such as rib waveguides, 



^ ■ can display nonlinear polarization effects in which the nonlinear phase difference dominates the linear 

O 

contribution, provided the birefringence is sufficiently small. We demonstrate that self-flipping polariza- 



tion states can appear in such rib waveguides at low (mW) power levels. We describe an optical power 



c/3 , limiting device with optimized rib waveguide parameters that can operate at low powers with switching 

O 
•'-^ ■ properties. 

I. Introduction 

►^ I Nonlinear interactions between the two polarization modes of a waveguide lead to intriguing 

physical effects and opportunities for new devices for optical data processing [|I]-[|4l. We show 



\^ I here how nonlinear interactions of this form can be employed to develop all-optical devices at 



low power levels, using switching (bistable) properties of the polarization phase difference, as 
described in [5|, at zero birefringence wavelength in subwavelength rib waveguides. In general, 
the interaction between the two polarizations in an optical waveguide is described by the coupled 
rS ■ Schrodinger equations 0): 



n=l 

= i (7,- \Ajf + 7, \Akf) A, + h',A*Al exp{-2izAPjk), 
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where j,k = 1,2 (j 7^ k) index the two polarization modes, Ai, A2 are the amplitudes 
of the corresponding fields, /3j„ are the nth order propagation constants, A(3jk = — A/3fcj is 
the linear birefringence, 7j,7c and 7^ are the effective nonlinear coefficients representing self 
phase modulation, cross phase modulation and coherent coupling of the two polarization modes, 
respectively. 

The weak guidance approximation assumes that the effective mode areas of the two polar- 
ization modes are equal [4], leading to 71 = 72 = 37c/2 = 87^. We have shown in [[51, |l6l 
that these equalities are not necessarily valid for waveguides with large index contrast and 
subwavelength dimensions. For such cases a new class of time-independent polarization states 
appear, in which the phase difference between the two polarizations oscillates abruptly between 
two well-defined values, which we refer to as switching (bistable) behavior [5]. Associated with 
this, the polarization state of the propagating mode flips abruptly through fixed angles. This is the 
result of competition between the linear and nonlinear phase differences of the two polarization 
modes propagating along the waveguide. While the linear phase difference is proportional to 
zAf3 (where (3 = (3i2 = —/32i), the nonlinear phase difference at any point along the waveguide 
depends on the 7 values and on the power coupled into each polarization. As a result, one 
expects that significant nonlinear polarization effects should be observed as A/3 approaches 
zero, i.e. when the linear phase difference is negligible. In waveguides with highly symmetric 
geometries, however, we have A/3 = only if there is also a symmetry between the polarization 
mode distributions (for example, in waveguides with circular or square cross sections), in which 
case the nonlinear phase difference also approaches zero for the same power coupled into the 
two polarizations. For this reason, we previously found polarization switching only for powers 
at the kW level [Sj. 

II. Rib waveguides 

Here, we consider rib waveguides with structural dimensions chosen such that the birefringence 
is nearly zero at the operating wavelength of 1550nm, following [|3-[|9l. For such cases, the 
waveguides are polarization independent i.e., A(3 = 0. Previous studies have identified the 
relation between the width and the height (etch depth) of the rib for independent polarization 
operation of rib waveguides [|71, [|9l and have used them in applications such as ring resonators 
[[HI. For these waveguides, the mode field distributions for the two polarizations are not identical. 





Fig. 1. Mode profiles for a rib waveguide with dimensions as shown. 



which leads to 71 7^ 72- Therefore, the nonlinear phase difference is large relative to the linear 
contribution, and leads to switching behavior at power levels several orders of magnitude lower 
than before [5J. Hence, we are able to demonstrate polarization self-flipping through 7r/2 at 
just mW powers. This approach produces polarization states which flip under their own power, 
whereas in [j2l, [|3l a second pump is required. 

An example of mode field distributions for two polarizations with A/3 = is shown in 
Fig. [H for a chalcogenide glass-based rib waveguide. Evidently the distributions for the two 
polarizations are not identical, which leads to 71 7^ 72. These profiles have been calculated using 
a finite element package, and the 7 values have been determined using the formulas (3) in [5]. 
The birefringence can be calculated as a function of the rib width and rib height (i.e. etch depth), 
and is plotted in Fig. [21 where the white line indicates points for which A/3 = 0. 

We now investigate nonlinear polarization as described in [l5]|, |l6l, for waveguides with A/3 ^ 
(near the white line in Fig. |2l left). We substitute Aj = ^JP'^e^'^i for j = 1,2 into Eqs. ([U), 
where Pj is the power of the field Aj with phase 0^ and look for time-independent solutions. 
The total power Pq = -Pi + -P2 is constant in z. We define the following dimensionless variables: 



-To 

A/3 



= 2A0, 
7c -72 



r = 2y^PoZ, 

I. 71 + 72 



27c 



(2) 
(3) 



where A0 = 4>i — 4>2 + zA/3 is the phase difference between the two fields. From Eqs. (d) 
we obtain a system of equations for v, 9 (see Eqs. (5,6) in [O) which can be solved for any 
given vq = v{0),9o = 9{0) at the initial location r = (with < i^o < !)» either analytically. 



or numerically for any specified values of a, b. There are three types of solutions: steady state 
solutions, periodic (including switching) solutions, and soliton solutions. For periodic solutions 
the dimensionless period T^. of v, cosO (regarded as functions of r) depends only on a, b and on 
Vq, Oq. Hence, the functions Pi, P2, cos2A0 are periodic in z with a period T^ = T^- / (2'j'^Pq) , as 
follows from the definition of r in Eq. ^. The switching behavior discussed in flU, ^ requires 
1 < a < 2b — 1 which from the definitions (|3]) implies that 



Poilc + 7c - 71) < A/3 < Po(72 - 7c - 7c)- 



(4) 



This inequality is satisfied at relatively low powers Pq provided A/3 is sufficiently small. 

Let L be the length of the waveguide, hence ^ z ^ L. We are interested in exploring the 
possible values of v, 9 at the endpoint z = L, i.e. at r = 2'y'^PoL; v, 9 each depend nontrivially on 
the total power Pq, regarded as a variable parameter, through a. One can in principle determine 
the explicit dependence of v, 9 on Pq with the help of the exact solutions, but for waveguides 
with A/3 = 0, however, both parameters a, 6 are fixed and in this case Pq enters the defining 
dimensioned equations as a scale factor only. This means that the periodic functions t>,cos6' 
and their period T^. are determined entirely by the values of the 7 coefficients, together with the 
initial values fo,^o- The functions Pi(z) , P2{z) , cos A(f)(z) at z = L can therefore be stretched 
or compressed by varying either Pq or L, in particular cosA(/) at z = L may be regarded as 
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Fig. 2. Contour plot of A/3 in m ^ (left), and the efficiency factor t] (right) defined in Eq. ^, for a range of rib waveguide 
parameters. 



a function of either L or Pq and, since 2A0 = 0{2'y'^PQL), is periodic in both L and Pq- The 
corresponding periods T^ and Tp^ of cos 2 A0, of dimension length and power respectively, are 
given by 

Tl = T, = 7^, Tp, = ^, (5) 



27^Po' ^" 27^L' 

where the dimensionless period T^ is independent of both L, Pq. These simple dependencies, a 

consequence of the scaling properties, are useful for the application outlined below. 

We consider now light that is initially linearly polarized, i.e. 9q = 0, with equal power coupled 
into the two polarizations, hence vq = 1/2. As an example, as explained further below, we choose 
a rib width of 990 nm and a rib height (i.e. etch depth) of 900 nm, with a substrate width of 
2400 nm and a height of 340 nm. We find, in units of (W.m)-^: 71 = 81.6,72 = 82.2, 7^ = 
49.6,7^ = 24.7, with numerical values for A/3 such that \A/3\ ^ 3.8 x 10~^ m^^ With 1 mW 
input power and aim long waveguide, the nonlinear phase change is more than 200 times 
larger than the linear contribution. 

With these 7 values we have a = 1.32,6 = 1.308, and so 1 < a < 26 — 1 is satisfied, and 
therefore switching states exist. Fig. |3] shows Pi,P2,cosA0 as functions of 2; at Pq = 1 W. 
Evidently Pi + P2 is constant, and Pi, P2 each have the period T^ = 0.47 m; cos A0 displays 
switching by flipping between ±1 with a period 2Tz = 0.94 m. 
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Fig. 3. Pi (blue), P2 (green) on the left, and cos A<j!> (blue), \dz cos A(/)| (red) as functions of z. 



We optimize the rib waveguide parameters by choosing values that result not only in A/3 = 0, 
but also such that cosA(/) flips as abruptly as possible between ±1, relative to the period. In 



order to quantify this transition, we define a transition length Ltrans over which cosA0 flips 
between ±1, and a dimensionless efficiency factor r] by 

2 T. 



^ trans 



V = -r^. (6) 



where 2Tz is the period of cosA0(z). The derivative d^ cos ^(l){z) is evaluated at cosA0 = 
since this is where the maximum slope of cos A0 occurs. This is evident in Fig. |3] (right), which 
plots |(92 cos A0(2;)| (red) as a function of z. 

For low power operations and fast polarization transitions we minimize Lt^ans relative to the 
period 2Tz, i.e. we maximize r] within the physical limit of a unit length waveguide, by varying 
the rib width and height but keeping A/3 = 0. This is equivalent to minimizing \a — 6|, which 
ensures that the initial value Vq = 1/2 lies close to the unstable steady state at f = (a— 1)/2(6— 1), 
as in the example shown in Figs. 5(ii) and 6(ii) of |[6l. Fig. [2] (right) shows 7] for a range of rib 
waveguides of varying rib width and rib height, where A/3 = holds on the white line. The 
maximum value of i] occurs at the top right of the white line, with larger values corresponding 
to periods 2Tz that exceed the length L of the waveguide. Choosing L ?^ 1 m, we arrive at the 
rib width 990 nm and rib height 900 nm as considered above, with r] = 12.3. 

in. Power limiting optical device 

The switching properties of cos A</), considered as a function of either Pq or the waveguide 
length L, enable us to construct in principal a power limiting optical device (a "surge protector"). 
Consider a rib waveguide with the structural parameters of the given example, with a polarizer 
at the output (at z = L) which is aligned either parallel or perpendicular to the direction of 
polarization of the input beam, which is linearly polarized. Denote the vector field at z = L hy 
Aout = ^ix + A2y, where x, y are the unit vectors along the X, Y axes respectively, and let 
r = (x + y)/v^ be the unit vector defining the direction of the polarization axis of the polarizer. 



The polarizer therefore projects the output vector field onto Aout ■ r = ^/Pi/2e'^^'^ + A/i\72 6*"^^ 
with a power output Pout = I Aout ■ rp- As a function of the initial total power Pq we have 



Pout(-Po) = -Po/2 + Pq V ^(1 ■" v) cos{6 /2) where f , 6 are known functions of r = 27^PoL. 

We choose L = Li = 0.4 m, so that switching behavior appears at power levels Pq ~ 1 W. 
The output power P^^^^ is shown in Fig. |4] (blue) as a function of Pq, and depends linearly on Pq 
for Pq ^ 0.3 W, but falls to zero in the range 0.9 — 1.4 W. Hence the polarizer in conjunction 



with the switching properties reduces the power output to zero over this range, and yet has no 
effect for small powers, i.e. P^^^{Pq) = Pq for Pq ^ 0.3 W. 




0.5 1 1.5 

Input power (W) 



CD 

o 



^0 



O 




0.05 0.1 0.15 0.2 
Input power (W) 





0.5 1 1.5 

Input power (W) 



Fig. 4. The power output Pout for waveguide lengths 0.4 m (blue) and 0.8 m (green), each with polarizers at the output end, 
and the combined two-stage output (red). 

In order to extend the range over which the output power is limited, we consider a second 
waveguide with the same structural parameters as before, but of length L2 = 2Li = 0.8 m, 
with an output P^^^ as shown in Fig. |4] (green); this plot is a compressed copy of the plot for 
Pout^ as follows from the scaling properties with respect to Pq and L. This second waveguide 
and polarizer are attached to the output end of the first waveguide, and the first polarizer then 
ensures that the initial conditions for the second waveguide are Oq = 0,fo = 1/2. The total 
output of the device as a function of Pq is given by -Pout(-Pout(-Po))5 which is shown in Fig. |4] 
(red). The output is linear for Pq ^ 0.15 W, but is limited to values ^ 0.2 W for all input powers 
Pq ^ 1.7 W. The maximum power operation range can be adjusted by varying the lengths and 
hence periods of the connected waveguides, and also by adding further stages. This conceptual 
design does not require any additional pumps for its operation, which is characteristic of the 
structurally-induced nonlinear polarization effects discussed previously [|3, [|6l. 



IV. Conclusion 

We have demonstrated the existence of nonlinear phase effects for waveguides with reduced 
symmetries, which are much larger than linear effects, at mW power levels. We have also 



demonstrated that switching behavior with 7r/2 polarization flipping can occur at the 1 W 
input power level. Although we have analyzed properties of the system for the case of zero 
birefringence, this is convenient but not essential; the exact solutions described in [|5l, [l6l allow 
one to determine the precise dependence of cos A0 on the input power for any A/3, including 
small nonzero values. The vectorial nonlinear model for sub wavelength waveguides has enabled 
us to uncover these phenomena. Our results show that applications displaying this nonlinear 
polarization behavior can, in principle, operate at practical power levels. 
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